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Response of an interferometer becomes complicated for gravitational wave shorter than the 
arm-length of the detector, as nature of wave appears strongly. We have studied how parameter 
estimation for merging massive black hole binaries are affected by this complicated effect in the case 
of LISA. It is shown that three dimensional positions of some binaries might be determined much 
better than the past estimations that use the long wave approximation. For equal mass binaries 
this improvement is most prominent at ~ 10 5 Mq. 



I. INTRODUCTION 

The Laser Interferometer Space Antenna (LISA), a joint project of NASA and ESA, is planed to be launched 
in 2011 B). LISA is sensitive to low frequency gravitational waves at 10 _4 Hz < / < 10 _1 Hz, and would directly 
confirm gravitational radiation from some known Galactic binaries Thousands of closed white dwarf binaries 
would be also detected by LISA. They can be regarded as insurances of the project, but other exciting phenomena 
might be observed 0. Coalescence of a massive black hole (MBH) binary is the most spectacular and energetic event. 
Gravitational wave from a MBH binary is in itself very interesting for studies of general relativity, but would also 
bring significant impacts on astrophysics and cosmology |2| , even though event rate is highly unknown at present [0 . 
Parameter estimation errors are basic measures for discussion of gravitational wave astronomy. For example, with 
small error box for the position of a MBH binary we might specify its host galaxy and investigate it with various 
observational tools. 

LISA is constituted by three space crafts that keep a triangle configuration with side (arm-length) I = 5.0 x 10 6 km 
and trail ~ 20° behind the Earth ||. Response of a detector becomes complicated for gravitational waves with 
wave-length A < I, as the nature of wave appears strongly This interesting feature might be an advantage 

for signal analysis but has not been taken in properly so far. In this paper we discuss how the parameter estimation 
errors are changed by this effect. 

This paper is organized as follows. In Section II we study the effect of a finite arm-length on the phase shift 5(f) 
of the detector, and introduce the past analysis that is based on the simple picture using the variations of the arm- 
lengths 51 under the long wave approximation. In Sec III A gravitational wave forms from MBH binaries are briefly 
discussed. Noise spectrum for both the variations 51 and the phase sifts 5(f) are presented in Sec III B. The spectrum 
for the latter is simply reproduced from that for the former. In Sec IV our numerical results are shown, and signal to 
noise ratio and the parameter estimation errors are compared for the two methods (with and without the long wave 
approximation). Though we mainly study gravitational waves from MBHs, results for nearly monochromatic sources 
(from e.g. Galactic binaries) are also presented. Sec V is devoted to discussion. In Appendix A we make same kind 
of comparisons for the noise canceling combination of data streams. In Appendix B we present explicit expressions 
for detector's response to gravitational wave from a binary. 



II. DETECTORS WITH FINITE ARM-LENGTH 

We name the three vertexes (space crafts) of LISA as A, B and C. Cutler studied parameter estimation of 
binaries from variations of arm-lengths in the following forms 

5l AB (t) - 5l AC (t), {(5l BA (t) - 5l BC {t)) - (5l CB (t) - 5l CA (t))} /V3. (1) 

The normalization factor l/\/3 in the second expression is explained later. Here 51 xy is the variation of the arm- 
length between two space crafts X and Y (X, Y £ {A, B, C}, X ^ Y). This simple picture for detector's response is 
valid only at long wave limit where gravitational wave-length is much larger than the detector's arm-length (see also 
HU). With quantities more close to detector's output these two data should be replaced by 

5<t>A B {i) - 5<f> AC (t), {(64> BA (t) - 5<f> BO {t)) - Wcb® - StcA®)} /V3, (2) 

where Sqbxy represents the single-arm phase shift of the laser beam that leaves a detector (vertex of the triangle) X 
for Y and then returns back to X. The quantity 5(f> AB — Sip A c is an interferometer signal obtained at the detector A. 
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In Ref . § it is assumed that the noise for the time variations SIxy is stationary, Gaussian and symmetric among 
three arms. We make same assumption for the corresponding phase shifts 54>xy In this case noises of the above two 
data (ph (or (El)) do not correlate to each other. 

Response of the phase shift 5cj>AB for gravitational wave h incoming from a direction fl s is expressed as follows 

' "' ) " Al:il > - 1 cos2i>AB[(l-cos6AB)h{t,A) + 2cosO A Bh(t-T,B)-(l + cosdAB)Ht-2T,A)}, (3) 



27w dt 2 



where Oab is the angle between the source direction fl s and the arm xa — xb, and vq is the fundamental frequency 
of the laser. tpAB is the principle polarization angle of the quadrupole gravitational wave h for the arm xa ~ xb- We 
denote the polarization basis tensor H a b of the wave as H a b — p a Pb — Qalb using two orthogonal vectors p and q with 
p ■ fl s = q ■ fl s = p ■ q = and |p| = \q\ = 1. Then the angle ipAB is given by t&nipAB = {q • {xa — xb))/(p • (xa — 
x b)) p0[ |- The quantity h(t,A) is value of the gravitational wave at point A and time t, and can be expressed as 
h(t + fl s ■ XA(t)/c) in the present case. We have denoted the propagation time of light for the arm I by r = l/c, and 
neglected very small relative motions between the vertexes in the time scale r. In Appendix B we present explicit 
expression of Eq.(^) for a binary source. 

When the gravitational wave-length is much larger than the arm-length I, namely / <C (2ttt)^ 1 = f arm , structure 
of wave is irrelevant at the spatial scale I and the wave h causes variations S<Pab or 81 ab in very simple manners. In 
this long wave limit the phase shift is directly proportional to the variation of the arm-length 81 xy as 84>xy /^ttvo — 
SIxy /2c, and two data sets (|l|) and (^) are equivalent. We have the following relation by perturbative evaluation of 
Eq.@ (see also Ref. go)) 

- snr Q A b cos2iph(t) = — — . (4) 



27T^oT 21 

LISA has arm-length I = 5.0 x 10 6 km, corresponding to the critical frequency farm — 0.01Hz. Note that the frequency 
/ does not appear in the above relation (Q). In this long wave limit the angular pattern function (usually denoted as 
F + ' x jnj) contains complete information of the angular dependence of the interferometer signal (e.g. 84>ab — 84>ac)- 
Two data streams (Q) are equivalent to responses of two 90°-interferometers rotated by 45° (see discussion in 0). 
The factor 1 / -s/3 in the second expression of ([l]) or (Q) is given for this purpose. 

When the gravitational wavelength c/ f is comparable to or shorter than the arm- length I, the phase shifts 8<pxY 
show complicated response that depends strongly on the frequency / and the source direction relative to the detector 
(see Appendix B). As discussed earlier, this troublesome but interesting effects have been cut down so far in analysis 
of parameter estimation errors for binary sources. Past analysis was performed in the following manner [R^^]. 

(i) The angular averaged transfer function T(f) between the phase shift 8<f> and wave amplitude h was used 
8(j> oc T(f)h to include amplitude modulation caused by finitencss of the arm-length I ^ 0. Its effect can be effectively 
absorbed in the noise curve ySh{f) for the wave amplitude h. 

(ii) Then the simple analysis with the variation of the arm-length SIxy (valid only at the long wave limit) was 
applied with the angular pattern function F + ' x for the angular dependence of the response (and also Doppler phase 
caused by velocity of detectors) . 

At high frequency limit / 3> farm the transfer function becomes T(f) oc / _1 due to cancellation of wave within 
arms. The measurement sensitivity of LISA for the phase shift 8<f> is \J Ss^(f) oc f° (in units of Hz -1 / 2 ) at 10 _2 Hz < 
/ < 10° Hz where the shot noise is dominant source of noise. Thus LISA has y/ Sh(f) oc y/~S/i$(f)T(f) 1 oc / at high 
frequency region [§;f|||. 

Here we discuss basic aspects of signal analysis with paying attention to the finiteness of the arm- length /. The 
most important quantity for detection of gravitational wave is the signal to noise ratio (SNR) (ll) . For its estimation 
the above mentioned method with the angular averaged transfer function would be effective for LISA considering its 
rotation and frequency average of chirping signals. But parameter estimation errors (e.g. error box for direction of a 
source) might become much smaller than the previous analysis when we include the complicated response of the phase 
shift as it is. Gravitational wave signal measured by a detector (as expression (|J)) is given by values of gravitational 
wave h at different positions and times. Time delay between them would also affect parameter estimation. These 
two effects, time delay and cancellation within the arms, are natural outcome of the finiteness of the arm-length and 
appear in a similar form, namely, linear combination of terms like exp[2Triffl s ■ (xa — xb)/c] in Fourier space of the 
time variable t. Our primary aim in this paper is to show how the parameter estimation would be changed if wc 
simply replace the arm-length variation SIxy (as in Eq.(|l])) with the corresponding phase shifts S(f>xY (as in Eq.(||)) 
without using the long wave approximation. 
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III. ANALYSIS OF GRAVITATIONAL WAVES FROM MERGING MASSIVE BLACK HOLE BINARIES 



A. Gravitational Waveform 

We investigate the parameter estimation errors expected in matched filtering method When the wave signal 
contains fitting parameters {7^}, their estimation errors (variances) for signal analysis are evaluated with using the 
Fisher's information matrix as (A-jiA-jj) = . We analyze the phase shifts (g) by extending earlier studies to 
include effects caused by the arm-length, but other points are almost identical to Ref. Q that uses the variations (|l|) 
under the long wave approximation. Detailed analysis for the phase shift is given in Appendix B. 

For gravitational wave emitted by a binary we adopt the restricted post-Newtonian approach, but higher-order 
harmonics could become important in some cases 113]. We use the 1.5PN phase for a circular orbit jjjj 



*(/) = 2irft c - <£ c - 
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where t c and <p c are integration constants, p is the spin-orbit coupling term, fj, z , M z and M. z are the reduced mass, 
the total mass and the chirp mass of the binary. All of the mass parameters are multiplied by the factor (1 + z) 

(z: redshift of the binary) with suffix z. The total mass M z is expressed by other two masses as M z — Ai\ \i z 3 ^ 2 . 

The post-Newtonian expansion parameter x is defined as x = {Gnc~ 3 M z f} 2 ^ 3 . We have used the stationary phase 
approximation for the Fourier transformed waveform. In our analysis total number of the fitting parameters is 10 as 
■M z , n z ,P, <pc t c , hi D (luminosity distance) and (6 S , (f) s ), (61, (f>i). The latter four parameters define the direction fl s 
and the orientation f2\ of the binary in a fixed frame at the barycentre of the solar system. In this paper we use the 
error ellipse 



AH 



i = 2tt sin 0i ^ (AO 2 ) (Atf) - (AO^) 2 , 



(6) 



for the angular parameters (i = s,l) |7j and AV = ADAfl s for the three dimensional position error of a source. We 
also fix f3 = for the true value of the parameter f3. 

Unless stated explicitely, we put the upper cut-off frequency f cu t of signal integration at fi SC o when the binary 
separation becomes r = 6M z /c 2 (l + 2), roughly corresponding to the inner most stable circular orbit. The observed 
wave frequency at this separation becomes 



fit 
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M, 
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Hz. 



(7) 



We start integration of the wave signal from the time when the binary is lyr before its coalescence in the observer's 
frame. The gravitational wave frequency /o at this starting time is given as follows 



fo = 1.9 x 10" 



0.87 x 10 5 Mf; 



-5/8 



Hz. 



(8) 



B. Noise Curve 



We make quantitative analysis using noise curve of LISA that is sum of the instrumental noise and the binary 
confusion noise. For the instrumental noise spectrum given for the wave amplitude h (almost equivalently for 51) we 
adopt the following function ]7j 

Si(f) = 5.049 x 10 5 M/) 2 + a 2 (f) 2 + a 3 (f) 2 } Hz" 1 , (9) 

where cti(f) = 10~ 22 ' 79 (//10 -3 Hz)~ 7 / 3 is mainly the acceleration noise, 012(f) = 10~ 23 - 04 is mainly the shot noise, 
and 013(f) = 10 _24,54 (//10^ 3 Hz) approximately represents the cancellation effects due to the finiteness of the arm- 
length. The last term reflects the angular averaged transfer function as explained before, and becomes important at 
/ > 10 _2 Hz (see e.g. Q). We also include the binary confusion noise S c (f) (estimated for one year observation) that 
is stronger than the instrumental one at / < 10~ 2,5 Hz and shows strong dependence on frequency / (l5). Around 
/ > 10~ 2 - 75 Hz the confusion noise decreases significantly as the Galactic binaries are resolved at higher frequencies. 
Its explicite expression is given as follows 
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, 10 -42.68 5/ -1.9 Hz -l ; / < 10 -3.15 j 

s c (f) = { lo- 60 - 325 / -7 ^- 1 , lO" 315 < / < 1(T 2 - 75 , (10) 
I io- 46 - 85 /- 2 - 6 ^- 1 , io- 2 - 75 < /, 

where / is written in units of Hz. 

Now let us simply recover the noise spectrum Ss,p(f) for the phase shift 54>ab — &4>ac from that for the variation 
51. At low frequency region / <C farm they are simply related to each other as expressions (Q) (||) and are essentially 
equivalent (see Eq.(Q)). We use the following functional shape for the instrumental noise of the phase shift Ss<f>(f) as 
Ss4>(f) oc (ai(/) 2 4- «2(/) 2 ) by removing the effect of the angular averaged transfer function from Si(f). Note that 
the confusion noise is expected to be much smaller than the instrumental one at the relevant frequencies / > farm- 
For the second data of expression (j|) we use the same noise curve as the first one. 

Our calculation using the phase shift is much more complicated than the past one. But our numerical results should 
coincide with these by past one when we decrease the arm-length I and use the shape of the noise curve Sh(f) instead 
of S$,p(f). We have confirmed that the parameter errors A 7; given in table 2 of Ref. (obtained by past simple 
calculation) is reproduced with this procedure. 



IV. RESULTS 



We have analyzed gravitational waves from MBH binaries with 300 realizations of random directions fl s and 
orientations 17;. Firstly, their SNRs and estimation errors of the fitting parameters {7^} are calculated in both (i) 
past approach (suffix 0) with angular averaged transfer function under the long wave approximation, and (ii) our 
method (suffix L) based on the phase shifts of expression (g). In figure 1 we present distribution for ratios of SNRs by 
two methods (SNR) L /(SNR) Q and position errors AV L /AV , (Af2 s ) L /(Af2 s ) Q and AD L /AD Q obtained for equal 
mass MBH binaries with redshifted masses M z = 10 4 + 10 4 , 10 5 + 10 5 and 10 6 + 10 6 Af Q . Note that these ratios do 
not depend on cosmological parameters or distance D to the MBH binaries. As shown in the bottom panel of figure 
1, ratios of SNRs are close to unity and difference of two methods is very small (within 20%). However, two methods 
show considerable difference for the position errors. For 10 5 + lO 5 A/0 MBH binaries the errors AVl become typically 
10 times smaller than the past estimate AVo. For some samples the errors Ay are reduced even by a factor ~ 10 -2 . 
Comparing errors Af2 s and AD, the former become smaller by the finiteness of the arm-length than the latter. But 
note that the angular resolution Af2 s is, roughly speaking, a product of two errors A9 S and A</> s as in Eq.(|^). 

Interestingly enough, difference between two methods is not a monotonic function with respect to BH mass and is 
most prominent around ~ 10 5 M©. This mass dependence can be understood as follows. At higher masses 3> 10 5 M Q 
the upper cut-off frequency fi SCO is smaller than the critical frequency farm- Therefore nothing is different between 
two methods with expressions ([!]) and With the quadrupole formula for gravitational radiation the time before 
coalescence is given by the frequency / and the chirp mass A4 Z as follows JTl[] 

/ f \ -8/3 / M \ ~ 5 / 3 

t G w = 8.4 x 10 4 4™ — sec. (11) 

G V 1 _2 ' 75 Hzy VO-87 x 1O 5 M 0/ / v ' 

For lower mass BH binaries (<C 10 5 M Q ) LISA moves longer than the arm-length I during the phase / > 10 _2 75 Hz 
where signal becomes very strong due to decrease of the binary confusion noise. Thus effective baseline of the detector 
becomes larger than the arm-length I — 5.0 x 10 6 km for smaller mass BH binaries and impact of the finite arm-length 
would decrease. 

Now we investigate various aspects of parameter estimation caused by the finiteness of the arm-length using a 
specific set of binary parameters. We pick up the binary that has the smallest volume ratio AVl/ AVq — 0.011 in 
300 realization of figure 1 for 10 5 + 10 5 M Q MBH binaries. It has angular parameters 8 S = 2.49, S = 0.03, di = 2.32 
and cj)i — 4.46. Here we present the estimation errors A7j, not ratios as in figure 1. To normalize the amplitude 
of the signal we take the redshift of MBH binaries at z = 1 with cosmological parameters fio = 0.3, Ao = 0.7 and 
Ho = 75km/sec/Mpc. In figure 2 estimation errors Af2 s , Afii, AD/D and A/i z /fi z are presented as functions of 
the upper cut-off frequency f cut that was fixed at fi SCO with Eq.(^|) in the case of figure 1. We fix the lower cut-off 
frequency at /o (eq.Q) 

We have found that the intrinsic binary parameters such as M. z ,ii z ,j3,t c and 4> c depend weakly on the cut-off 
frequency f cu t- However, errors for the binary position Ai7 s , AD/D or its orientation A17; decrease significantly at 
fcut ^ 0.01Hz where the long wave approximation breaks down. This shows remarkable contrast to the past analysis 
shown with thin lines. Our results seem reasonable, as the response of a detector with finite arm-length depends 
strongly on the direction of the source fl s , and information of the distance D or orientation 17; is tightly correlated to 
them (Appendix B). Significant reduction of position errors AVl at higher frequencies would give further motivation 
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for studies of nonlinear gravitational dynamics (e.g. Post-Newtonian approach) that is a very tough problem on 
general relativity. By analyzing gravitational wave close to the final coalescence we might identify the host galaxy of 
a MBH binary! 

Next let us make hypothetical experiments to clarify some interesting points. We use the same set of the parameters 
as in figure 2 with f cut given by Eq.(j^). As commented before, the finiteness of the arm-length causes two similar 
effects (i) cancellation of waves within the arm and (ii) time delay between the vertexes of LISA. To extract effects only 
of the latter we calculate the volume error AVl' using data from three interferometers that exist at three vertexes 
of LISA but have arm-length I — > with angular averaged transfer function for sensitivity of h. Thus only the 
positions (separation) of the detectors are different from the past analysis that take the separation I — 0. We obtain 
AVi'/AVo = 0.15. This result indicates that two effects work cooperatively. Next we stop motion of LISA and keep 
its position at time t = t c . In this case the volume error AVl» becomes AVz/'/AVo = 0.012 and is very close to 
AVl/AVq = 0.011 that includes motion of LISA. In the past analysis (with I = 0) we use the amplitude modulation 
through the pattern function and the Doppler phase modulation both caused by motion of LISA || , and cannot solve 
degeneracy of sources direction f2 s and other variables when LISA stops. The response of a detector with I ^ 
depends strongly both on angular variables and frequency of incoming waves at / -J farm (see Appendix B). Thus we 
can solve the degeneracy for chirping binaries even without motion of LISA, though there would be two solution for 
detector's signal due to the symmetry of source-detector configuration. This is a qualitatively interesting point. 

Considering the mass dependence of our results, it is expected that the past simple method using expression (Q) 
under the long wave approximation would be effective for studying galactic compact binaries. These binaries with 
/ < 0.1Hz are nearly monochromatic and have long durations tew 3* lyr m the LISA band due to their small chirp 
masses 0]. We have investigated binaries with mi = m 2 = 1-M©, and the time to coalescence tew = 30yr and lOOyr. 
The wave frequency becomes / = 0.071Hz at taw = 30yr and / = 0.045Hz at t GW = lOOyr. We fix the observation 
period at lyr, thus do not observe the final coalescence in contrast to the analysis for MBH binaries. We calculate 
parameter estimation errors in this situation with random direction fl s and orientation It is found that the 
differences between (AVt,AVo) or ((SNR)l,(SNR)q) are less than 15 percent. Recalling that we have made a simple 
treatment for the effects of the transfer function, this result seems excellent. 



V. DISCUSSION 



In this paper we have studied how data analysis of gravitational waves from MBH binaries are affected by finiteness 
of the arm- length of LISA with using the phase shift of detectors. We have numerically confirmed that the past method 
with the long wave approximation is very effective for estimation of SNRs that are the most important quantities 
for detection of gravitational waves. However, LISA is able to observe merging MBH binaries with significant SNRs 
(> 1000) even at cosmological distances ^,0^,^] , and aspects of gravitational wave astronomy are more relevant for 
them, rather than SNRs. We have examined the parameter estimation errors Aji expected in both two methods, and 
shown that three dimensional position of MBH binaries could be determined much better than the past estimations. 
In the case of equal mass binaries differences between two methods are most prominent at ~ 1O 5 M and the volume 
of the error box can decrease significantly. We have shown that (i) this reduction is mainly cause by gravitational 
waves close to the final coalescence and (ii) position of chirping binaries can be in principle, estimated even without 
motion of LISA. The former would give further meaning for studies on strong gravitational dynamics, and the latter 
makes remarkable contrast to former discussions. 

The author would like to thank an anonymous referee for helpful comments to improve the manuscript. This work 
was supported in part by Grant-in- Aid of Scientific Research of the Ministry of Education, Culture, Sports, Science 
and Technology No. 0001416. 



APPENDIX A: NOISE CANCELING COMBINATION 

In the main text we have basically followed the model for signal analysis in Ref. JtJ , and simply replace the variations 
(|l|) with the corresponding phase shifts (|^) without resorting to the long wave approximation. It has been discussed 
recently that the laser phase noise can be removed well by devising combination of data at different times [fl6| . This 
noise is caused by un-equal arm-lengths of space detectors in contrast to the ground-based ones. In this appendix 
we discuss parameter estimation using this combination. Let us consider the situation that the space crafts B and C 
coherently transmit laser beams back to the space craft A. Then the following data stream X(t) is a noise canceling 
combination 

X(t) = 54> AB {t) - 84> AC {t) - 54> AB {t - 2l AC /c) + 5^ AC {t - 2l AB /c). (Al) 
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For quantitative study we take the limit Iab — I AC — I — 5.0 x 10 6 km. Then the Fourier transformation of X(t) is 
related to that for 5<pAB (*) — 8(j>Ac(t) as 



X(f) = (1 - eM^itthWMf) - HacU))- (A2) 

In this case the noise curve for the gravitational wave amplitude with data X is identical to that with &4>ab — 34>ac 
H , if we only include the acceleration and shot noises. 
In figure 3 we present the two ratios (SNR)l/(SNR)q and AVl/AVq for the noise canceling combination X given 



in Eq.(Al). Note that the factor (1 — exp (2-zr il f / c)) in Eq. ( |A2| ) contains none of our ten fitting parameters. Thus 
results for the data X(t) is identical to that for the single data stream 54>ab — &4>ac- Figure 3 shows same kind of 
mass dependence as figure 1. Difference between AVl and AVb is most prominent at mass ~ 10 5 M Q again. But the 
ratio AVl/ AVb is generally close to unity and effects of the finite arm-length are smaller. 

Note added After submitting this paper, there appears [jl7| that discusses three data streams A, E, T whose noises 
do not correlate. Our study can be easily extend to these data and more realistic results would be obtained. 



APPENDIX B: PHASE SHIFT OF A DETECTOR FOR GRAVITATIONAL WAVE FROM BINARIES 

In this appendix we give an explicit expression for the phase shift in the form (eq. [3] ) 

1 dS(f> AB {t) 



dt 



(Bl) 



The shift 8(f) ab is defined for the laser beam that leaves the detector A for B, and then returns back to A. For source 
of gravitational radiation h we consider a chirping binary at luminosity distance D, direction fl s , and orientation Q\. 
We denote xy as the position vector of a detector Y £ {A, B, C}, and use a coordinate system fixed to the barycentre 
of the solar system. As we do not directly use the angular pattern functions F + and F x (see fn}), there is no need 
to introduce a coordinate system fixed to detectors. 
The inclination angle i of the binary is given as 

cosi = fl s ■ ili. (B2) 

We denote the propagation time of light for the arm-length by r = \xa — xb\/c, and define the angle 9ab between 
the direction of a binary fl s and the arm xa — xg as 

fl s ■ (x A - x B ) 

cos Bab = ■ (B3) 

CT 

It is convenient to use the principle axes (p, q) for analysis of gravitational wave from a binary. These two vectors 
are expressed in terms of direction fl s and orientation fli as 

With these vectors gravitational wave from a binary is decomposed to two polarization (+ and x ) modes whose phases 
differs by n/2. The plus (+) mode has polarization basis tensor H^ b = p a Pb — QaQb and the amplitude A + at the 
Newtonian order as 



or 5 / 3 A/f 5 / 3 

A + = G n , (^) 2 / 3 (l + cos 2 »), (B5) 
Dc 1 



where M. z is the redshifted chirp mass of the binary. For this mode the principle polarization angle tpAB is given by 



, , {x A -x B )-q m „, 

tan^s = ^ . (B6) 

(x A - x B ) p 



The cross ( x ) mode has the amplitude 



x = 7 cos*. ( B7 ) 



G 



with the polarization tensor iff. = p a qb + g a Pb . Its principle polarization angle differs by 7r/4 from the plus mode. 

Now we can write down the phase shift ( Bl ) with various parameters of the binary. For notational simplicities we 
define a function U (t, Y) that contains information of the phase of the gravitational wave at the detector Y and time 
t (taking its origin appropriately) as 

U(t, Y) = cxp[2TTif(t + x Y ■ O s /c)\. (B8) 
Then the quantity V is the real part of the following expression (Eq.(3)) 

- (cos 2ip AB A + + i sin 2ip AB A x ) 
x [(1 - cos 6 AB )U{t, A) + 2 cos 9 AB U{t — t,B) — (1 + cos 9 AB )U(t - 2r, A)] . (B9) 

The expression in the square bracket of the above result can be written as 

u(t,A)R(6 AB Jr), (BIO) 

where we have defined the factor R as 

R{9 ab Jt) = [(1 - cos 9 AB ) + 2 cos 9 AB exp {2mfr(- cos 9 AB -1)}-{1 + cos 9 AB ) exp{-4mfr)] . (Bll) 

It is a simple task to obtain the Fourier transformation of the phase shift 8(f) ab- In a similar manner we can make 
the expression for the phase shift of another arm e.g. S(j) A c. Then we obtain the signals such as 5<fi AB — 54> A c- 

When gravitational wave-length is much smaller than the arm-length (/r ^S> 1), the factor R shows complicated 
response that depends strongly on both the angle fl s and the frequency /. In the long wave approximation (/r <C 1) 
we can perturbatively expand Eq.(Bll), and obtain 



R = Amfrsm 2 9 AB . 



(B12) 



Now the factor R becomes very simple. 

Parameters such as A4 Z , [i z , P are closely related to the time evolution of the frequency (chirp signal) due to 
gravitational radiation reaction as given in Eq. (^) . Thus the complicated response of the factor R has weaker impact 
on estimation of these parameters than the direction of the source f2 s (see figure 2). Distance D or the orientation 
fli of the binary a re d e termined by usin g the information of the amplitudes cos2yb AB A + and sin 2ip AB A x . As seen 
in the expressions (B2)(|B4|)(B5)(B())(B7), they are strongly correlated to the direction J7 S of the binary. Thus their 
estimation errors are also improved by the effect of the factor R. 
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FIG. 1. Distribution of relative magnitude of the three dimensional error box AV, SNR, the angular resolution Af2 s , and 
the distance error AD/D. We compare results from two data streams (with suffix L) and those obtained by past simple 
method (with suffix 0). The solid lines are results for MBH binaries with redshifted masses 10 6 + 1O 6 M , the doted lines for 
10 5 + 1O 5 M , and dashed lines for 10 4 + 1O 4 M . We have analyzed 300 binaries with random directions and orientations. 
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FIG. 2. Dependence of the parameter estimation errors A 7; on the upper cut-off frequency f cu t < fisco- We star integration 
of the signal from fo = 1.9 x 10~ 4 Hz. The MBH binary has redshifted masses 10 5 + 1O 5 M0, and exists at z = 1 with 
direction 9 S = 2.49, <f> s = 0.03 and orientation 8 3 — 2.32, <j) a — 4.46. The thin lines are the past estimations and thick ones 
are the new estimations. The solid lines represent for AD/D, the long dashed lines for Afi;, the short-dashed lines for A£l s , 
and the dash-dotted lines for the reduced mass Afi z /fi z . SNR becomes ~ 1054 for f cut = 7.0 x 10" 3 Hz and ~ 1140 for 
f cut = 2.0 x 10- 2 Hz. 
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FIG. 3. Same as figure 1, but for the noise canceling combination X. 
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